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^ : Abstract 

~r\ ' A relativistic mean field description of collective excitations of atomic nuclei 

is studied in the framework of a fully self-consistent relativistic random phase 
approximation (RRPA). In particular, results of RRPA calculations of multi- 
pole giant resonances and of low-lying collective states in spherical nuclei are 
analyzed. By using effective Lagrangians which, in the mean-field approxima- 
te ' tion, provide an accurate description of ground-state properties, an excellent 
agreement with experimental data is also found for the excitation energies of 
low-lying collective states and of giant resonances. Two points are essential 
for the successful application of the RRPA in the description of dynamical 
properties of finite nuclei: (i) the use of effective Lagrangians with non-linear 
terms in the meson sector, and (ii) the fully consistent treatment of the Dirac 
sea of negative energy states. 
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I. INTRODUCTION 

Relativistic Mean Field (RMF) models based on effective Lagrangians with non-linear 
terms in the meson sector |1]||, have been successfully applied in the description of static 
properties of nuclei. With only six or seven phenomenological parameters (coupling con- 
stants and masses), carefully adjusted to a small set of experimental data in spherical nuclei, 
excellent results are obtained for binding energies and radii, density distributions, deforma- 
tion parameters and moments of inertia over the entire periodic table ||. There are, of 
course, cases where the mean field description breaks down, and certainly there will be fur- 
ther refinements of the parameter sets in the future, but in general the relativistic mean field 
framework provides us with an unexpectedly simple and successful tool to describe static 
properties of the nuclear many-body problem. 

The description of excited states within the same framework has been much less studied. 
Very successful applications of the cranking relativistic mean field model in the description 
of collective rotational bands have been reported j|,|5| , and an excellent description of exci- 
tation energies of multipole giant resonances has been obtained with the Time-Dependent 
Relativistic Mean Field (TDRMF) model H0]. The cranking RMF model is, however, quasi- 
static, and in some cases the large computational effort prevents an accurate description of 
excited states in the time-dependent framework. In addition, with the exception of the 
monopole case, time- dependent calculations violate rotational symmetry, and therefore it is 
not always easy to separate modes with different angular momenta. 

In the non-relativistic mean field framework the standard tool for the description of 
excited states is the Random Phase Approximation (RPA) [[sHlOU- On the other hand, 
applications of the Relativistic RPA (RRPA) have only very recently reached quantitative 
agreement with experimental data, and so far only in a few very specific cases |IT| , [12| . There 
are two reasons. First, although RRPA computer codes have been developed already in the 
eighties |T3Hl7|, those were based on the linear version of the relativistic mean field model, 



sometimes also called the Walecka model. With only linear terms in the meson sector, the 
RMF model does not reproduce ground state properties of finite nuclei on a quantitative 
level, and therefore also the corresponding RRPA calculations of excited states are not 
in agreement with experiments. Second, only recently it has been shown that the fully 
consistent inclusion of the Dirac sea of negative energy states in the RRPA is essential for a 
quantitative comparison with experimental excitation energies of isoscalar giant resonances 



||19|| . It has to be emphasized that already in Ref. |L8| it was shown that the inclusion of 
configurations built from occupied positive-energy states and empty negative-energy states 
is essential for current conservation and the decoupling of the spurious state. What was not 
known until very recently, however, is that the inclusion of states with unperturbed energies 
more than 1.2 GeV below the Fermi energy has a pronounced effect on giant resonances 
with excitation energies in the MeV region. 

In the last decade new techniques have been developed which enable the inclusion of 
non- linear meson interactions in the RRPA |2D| . With the use of non-linear Lagrangians it 
became evident that the a/i-configurations (a empty state in the Dirac sea, h occupied state 
in the Fermi sea) strongly affect the excitation energies of the isoscalar giant resonances. 
The a/i-configurations have recently been included in the calculation of the relativistic linear 



response [11 1, and in the diagonalization of the RRPA-matrix [12]. The pronounced effect 



of the Dirac sea in the RRPA calculation of giant resonances was completely unexpected, 
because it occurs in particular for states (T = 0, J n = + or 2 + ) which are orthogonal to 
the spurious mode (T = 0, J n = 1~). It is caused by the fact that, due to the relativistic 
structure of the RPA equations, the matrix elements of the time-like component of the vector 
meson fields which couple the a/i-configurations with the p/i-configurations, are strongly 
reduced with respect to the corresponding matrix elements of the isoscalar scalar meson 
field. As a result, the well known cancellation between the contributions of the a and uj 
fields, which, for instance, occurs for the ground-state solution, does not take place and 
large matrix elements which couple the a/z-sector with the p/t-configurations are present in 
the RRPA matrix. In addition, the number of a/t-configurations which can couple to the 
p/i-configurations in the neighborhood of the Fermi surface is much larger than the number 
of p/i-configurations. Because of level repulsion (Wigner's no-crossing rule) the effect is 
always repulsive, i.e. the inclusion of the a/i-configurations pushes the resonances to higher 
excitation energies. 

The non-linear meson terms in the effective Lagrangian, and the configurations built from 
negative energy Dirac sea states, were included in the time-dependent calculations of Refs. 
IK?]. A rigorous derivation of the RRPA-equations in the small amplitude limit of TDRMF 
PI , shows that the time-dependent variation of the density matrix 5p contains not only the 
usual p/i-components with a particle state p above the Fermi sea and a hole state h in the 
Fermi sea, but also components with a particle a in the Dirac sea and a hole h! in the Fermi 
sea. Although for the stationary solutions the negative-energy states do not contribute to the 
densities in the no-sea approximation, their contribution is implicitly included in the time- 
dependent calculation, even if the time- dependent densities and currents are calculated in the 
no-sea approximation. The coupled system of RMF equations describes the time-evolution 
of A nucleons in the effective mean-field potential. Only particles in the Fermi sea are 
explicitly propagated in time, and the corresponding vacuum is, of course, time-dependent. 
This means that at each time there is a local Fermi sea of A time-dependent spinors which, 
of course, contain components of negative-energy solutions of the stationary Dirac equation. 
One could also start with the infinitely many negative energy solutions and propagate them 
in time with the same Dirac hamiltonian. Since the time-evolution operator is unitary, the 
states which form the local Dirac sea are orthogonal to the local Fermi sea at each time. 
This is the meaning of the no-sea approximation in the time- dependent problem. 

Both effects, the non-linear meson terms in the effective Lagrangian, and the inclusion of 
Dirac sea states in the RRPA, are now well understood. In this work we present a system- 
atic study of different excitation modes in the fully consistent RRPA framework. By using 
effective Lagrangians which, in the mean-field approximation, provide an accurate descrip- 
tion of ground-state properties, we will show that an excellent agreement with experimental 
data is also obtained for the excitation energies of low-lying collective states and of giant 
resonances. 

The paper is organized as follows: the formalism of the fully consistent RRPA is presented 
in the section II. In section III we analyze in detail the lowest multipole isoscalar and isovector 
excitation modes, both at low energy and in the regions of giant resonances. The results are 
summarized in section IV. 



II. THE RELATIVISTIC RANDOM PHASE APPROXIMATION 

In the RRPA framework two methods are used to calculate the physical quantities rele- 
vant to nuclear excitations. Either the linear response function of the nucleus is calculated 
directly, or the RRPA matrix equation is diagonalized. In both cases the starting point is 
the self-consistent solution for the ground state. The Dirac equation and the equations for 
the meson fields are solved by expanding the nucleon spinors and the meson fields in terms of 
the eigenf unctions a spherical harmonic oscillator potential 0] . The resulting single nucleon 
spinors are used to build the RRPA basis: a set of particle-hole (ph) and antiparticle-hole 
(ah) pairs. The number of basis states is determined by two cut-off factors: the maximal 
p/i-energy (e p — €h < E max ), the minimal a/i-energy (e a — e^, > E min ), and by the angular 
momentum coupling rules. In this basis the unperturbed polarization operator 11°, or the 
RPA matrix is built. Both methods have been employed in the present analysis, and it has 
been verified that they produce identical results. 

The linear response of a nucleus to an external one-body field is given by the imaginary 
part of the retarded polarization operator 

R(F, F; k, w) = -ImU R (F, F; k, k; u) , (1) 

7T 

where F is the external field 4x4 matrix operator. The technical details of the calculation 
can be found, for instance, in Ref. [pOR . From the linearized Bethe-Salpeter equation in 
momentum space, the polarization operator is obtained 

U(P, Q- k, k', w) = U°(P, Q- k, k', w) + (2) 

J2qI f d 3 hd 3 k 2 u°(p, r 0; k, k 1; uOD^ki k 2 , cu)n(r , q- k 2 , k', u), 

where the summation is over the meson fields which mediate the effective nucleon-nucleon 
interaction. P, Q and T^ are 4x4 matrices which determine the Dirac structure of the 
corresponding operators and propagators. For the isoscalar scalar meson (0 = a) Y a = 1, 
for the isoscalar vector meson (0 = u) T = 7 for the time component, and T = 7 ck for 
the spatial components. An additional isospin matrix f has to be included for the isovector 
vector meson (<p = p). In general, the meson propagators are defined as solutions of 

^ 9 . + ^ri*^(^) = -S{x-y) } (3) 

where, for an effective Lagrangian with meson self-interactions, U^(4>) is the generalized 



mass term. For the non-linear a model [E6|,B7 



UM = \mla 2 + l -g 2 a" + \g,a\ (4) 



Effective Lagrangians with u self-interactions have also been considered ^8 



U u (u) = -m^^ + -c 3 (a;%) 2 . (5) 



The propagators in momentum space are obtained by the Fourier transformation of Eq. ([3D 
{to 2 - k 2 )Z^(k, k', w) - -^ j S^k - k^D^kt, k', to)d 3 h = (2vr) 3 5(k - k'), (6) 

where S<j,(k — ki) is the Fourier transform of gf 2 |</> > 

^(k-kO = /e-^^|^ 3 r. (7) 

Without meson self-interaction terms the propagators are diagonal in momentum space, 
D ^(kx k 2 , uj)— 5 (ki — k 2 )-D^(k 1 ) and the solution reads 

1 1 

Dj(k,u) = ±——-— — — r —— , (8) 

(2ir)- i uj 2 — k 2 — mi + %r\ 

where the plus (minus) sign corresponds to vector (scalar) mesons. Retardation effects (the 
frequency dependence of the propagator) are easily taken into account in the linear response 



formalism. This is not the case in the time-dependent formulation p2| , or for the solution of 
the RPA equations by matrix diagonalization in configuration space [|TI| . However, since the 
meson masses are very large as compared to excitation energies of nuclear collective modes, 
retardation effects can be neglected to a very good approximation. 
The unperturbed polarization is given by 

n°(r, r' ; k, k, u ) = y ^\ Ve ^ r \p)^\ v '^ r \ h ) _ (P\re-^\h)(h\T'^\ P ) 

~£ uj - e p + e h + ir] u + e p - e h + irj 

v (/i|re- ikr |a)(a|rV k ' r |/i) (a|r e - ikr |/i)(/i|rV k ' r |a) 



ah 



w - e a + e h + ir) u> + t a - t h + ii] 



where the indices p and h run over all unoccupied and occupied (Fermi sea) positive energy 
solutions of the Dirac equation, respectively. The index a denotes states in the Dirac sea 
(negative energy states). The single-particle matrix elements read 

(plTe^lh) = |d 3 r^t( r ) 7 o r e jkr ^(r) . (10) 

In spherical nuclei rotational invariance reduces the dimension of the basis. By using angular 
momentum coupling techniques, from the linearized Bethe-Salpeter equation (|J), separate 
equations are obtained for each angular momentum channel J. For a state with angular 
momentum J, the Hartree polarization reads 

n° ( r,r';M» =VV MPlj(k)\\f)(f\\P LIJ (kf)\\h) _ (f\\PlAk)\\h)(h\\P LIJ (V)\\f) ^ 

JXfX' uJ-e f + e h + ir] u + e f -e h + ii] 

(11) 

where the index / denotes particle states above the Fermi level (p), and negative energy 
states in the Dirac sea (a). Pljm and P' L i JM denote the spherical tensor operators 



P L jM(k) = j L (kr)[T x Y L {r)] JM and P' UJM {k) = jv{kr)[Y' x Y v (v)] JM . (12) 

The Dirac matrix V is coupled with the spherical harmonic Ylm to the total angular mo- 
mentum J . The reduced matrix elements are defined by 

^}(r) 7 ° jL (fcr)[r x Y L (r)] JM Mr) = R^ ( _^ ^ ^ ) (f\\pLj(k)\\h) . (13) 

For the scalar mesons and the time-like part of the vector meson fields, T carries angular 
momentum zero. In this case L = J. For the spatial part of the vector meson fields T = 7 a. 
The spin matrices u carry angular momentum 1, and therefore L = J — 1, J+l. The matrix 
elements of L = J vanish due to parity conservation. 



Alternatively, the RRPA matrices are built in configuration space [23.19 





A — ( ( C P ~ e h)S pp '5hh' | -I- | Vph'hp' Vph'ha' \ ,-,a\ 

\ (e a — eh)S aa t5hh> J \ V a h'hp' Vah'ha' J 

B = ( \f hh ' Xf hh ' ) (15) 

\ Vap'hh' Vaa'hh' J 

The two-body interaction V has the form 

V(r u r 2 ) = ^(n, r 2 ) (3^(3™ + D w (n, r 2 ) (l - cx^cx^) . (16) 

Here, ay' and (3^ are the Dirac matrices acting in the space of particle i, and for simplicity 
the contributions of the p-meson and the photon have been omitted. 
The eigenmodes of the system are determined by the RRPA equation 

(17) 

The transition densities and transition amplitudes are expressed in terms of the RRPA 
amplitudes (X U ,Y U ) p9[ . The dimension of the RRPA matrix is twice the number of ph- 
and a/i-pairs, and the matrix is obviously non- Hermit ian. In the non relativistic case the 
RPA matrix is real and the matrices (A ± B) are positive definite. The eigenvalue equation 
is then reduced to a Hermitian problem of half dimension (for details see Ref. [[^J). In 
the relativistic RPA, on the other hand, the matrix A has large negative diagonal elements 
(e a — e/i) < —1200 MeV, and the eigenvalue equation can no longer be reduced to a Hermitian 
problem of half dimension. The full non-Hermitian RRPA matrix has to be diagonalized, 
even in cases when the matrix elements are real. This, of course, presents a considerable 
increase in computation as compared to the non relativistic case. 

III. COLLECTIVE EXCITATIONS AND GIANT RESONANCES 

In the present analysis we have performed RRPA calculations with nonlinear effective 
Lagrangians which, in the mean-field approximation, describe ground-state properties of 
finite nuclei all over the periodic table. Within the non-relativistic framework, the RPA with 
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both the Hartree-Fock mean field and the residual particle-hole interaction derived from the 
same energy functional [|3HT0|1 has been very successfully applied in calculations of low-lying 
collective states and giant resonances. Some of the earliest application of the relativistic 
RPA have also used the residual particle-hole interaction derived from the the same effective 
Lagrangian which generates the solution for the ground state. These calculations [p~3| , p~6|] , 
however, were based on the most simple, linear a - u relativistic mean field model, and only 
a qualitative comparison with experimental data was possible. 



Most results in the present analysis have been obtained with the NL3 |M| non-linear 
parameter set of the effective Lagrangian. We have used this effective interaction in most of 
our recent applications of the relativistic mean-field model, and the results that have been 
obtained for the ground-state properties indicate that NL3 is probably the best effective 
force so far, both for nuclei at and away from the line of /3-stability. In order to study how 
the results for excited states depend on different effective interactions, we have also used 
other well known non-linear parameter sets: NL1 ||26||, NL-SH [Effl and TM1 P8|, as well as 



the linear Walecka model HS pqj . The different parameter sets are listed in Table 1. 

The multipole strength distributions have been calculated either with the response func- 
tion method or by diagonalizing the RRPA matrices. In the former case the strength function 
is given by Eq.O) in the limit of k — > 0. The isoscalar and isovector multipole operators are 
defined 

A 

F Xti = J2 r i Y \^iAi) for isoscalar 

A 

= y^ 73 r^Yx^Oj, 4>j) for isovector excitations. (18) 

4 = 1 

An energy averaging is performed by replacing the infinitesimal quantity irj — > iT, with 
T=2 MeV. In the matrix RRPA diagonalization method the resulting multipole transition 
probabilities between the ground and excited states 

B(n u ) = \(u\F\0)\\ 

are used to calculate the Lorentzian-averaged strength distribution 

*fi) = S*(",) (g _ffi +n/4 - 

r = 2 is the width of the Lorentzian distribution. 

A. Low-lying isoscalar states 

As an example of low-lying vibrational states in spherical nuclei, in Table 2 we display the 
RRPA energies and B(EL) values for the lowest collective states with angular momentum J n 
= 2 + , 3~, 4 + , 5~ in 208 Pb. The effective interaction is NL3 and the results have been obtained 
by diagonalizing the RRPA matrices in configuration space. The position of the low-lying 
collective states, as well as the B(EL) values, depend sensitively on the details of the single- 
particle spectrum around the Fermi surface. The general agreement with experimental data 
is very good, both for the excitation energies and transition probabilities. 



Similar results were obtained in the non-relativistic consistent RPA framework with 
Skyrme interactions ||, and with the Gogny interaction ||. 

B. Giant monopole resonances 

The study of isoscalar giant monopole resonances (ISGMR) in spherical nuclei is essential 
for obtaining information on the nuclear matter compression modulus K nm . Through the 
nuclear matter equation of state, this quantity defines basic properties of nuclei, supernovae 
explosions, neutron stars and heavy-ion collisions. Both the time-dependent relativistic 
mean- field model and the relativistic RPA | H| , have been recently applied in microscopic 



calculations of monopole excitation energies. The non-linear effective interactions have much 
lower nuclear matter compression moduli, as compare to older, linear parameter sets. It has 
been found that experimental GMR energies are best reproduced by an effective force with 
K nm « 250 — 270 MeV, i.e., about 10% - 15% higher than the values found in non-relativistic 
Hartree-Fock plus random phase approximation (RPA) calculations. The difference between 
the non-relativistic and relativistic predictions for the nuclear matter compression modulus 
is so far not understood, and it requires further investigations. 

There are less experimental data on the isovector giant monopole resonances (IVGMR). 
Because of the repulsive character of the particle-hole interaction in the isovector channel, 
IVGM resonances are found at higher excitation energies (typically 25 - 30 MeV in heavy 
nuclei). The widths of the IVGMR are also much larger (typically around 10 MeV) than 
those of the isoscalar resonances. In Fig.|l] we display the isovector monopole strength 
distribution in 208 Pb, calculated with the NL3 effective interaction. The large width of the 
calculated resonance is caused by Landau damping. Similar features were also found in non- 



relativistic RPA calculations with Skyrme interactions |3(J|. In Table 3 we have compared 
the centroid energies of the IVGMR in 208 Pb, calculated with different non-linear effective 
interactions, with experimental data. The agreement with experiment is very good, and only 
small differences are found between the predictions of the different parameter sets. This is, 
of course, expected since these effective interactions were adjusted to give similar values for 
the nuclear matter symmetry energy. 

C. Giant dipole resonances 

Except for the monopole mode, the response of a nucleus to higher multipoles is difficult 
to calculate in the time-dependent RMF. Rotational invariance is broken and the differential 
equations have to be explicitly solved at each step in time on a two dimensional mesh in 
coordinate space. It is difficult to keep the solutions stable for very long times, which are 
necessary for high accuracy. It is much easier to perform RRPA calculations for higher 
multipole modes. 

The isoscalar dipole mode has been recently studied in the RRPA framework HTl~l , |l"2f , and 
in the present work we analyze the well known isovector giant dipole resonance (IVGDR). 
The effect of the Dirac sea states on the isovector dipole strength distribution is illustrated 
in Fig.|^. The IVGDR response function in 208 Pb has been calculated with the NL3 effective 
interaction. The solid and dashed curves are the RRPA strengths with and without the 



inclusion of ah configurations, respectively. The effect of not including negative energy 
states is very small: the height of the peak is slightly lowered, and some extra strength 
appears in the 5-7 MeV region, but the energy distribution is practically unchanged. This 
is, of course, due to the fact that the ah configurations become really important only if the 



ex-meson dominates the RRPA response |Il|Jl9| . For the isovector modes the contribution of 
the isovector p meson is the dominant component of the particle-hole interaction. This is also 
illustrated in Fig.^, where the dashed-double dotted curve corresponds to the calculation 
with only the p meson field in the residual particle-hole interaction. Comparing the full 
RRPA response and the p-only calculation, it is seen that the only effect is that some 
strength from the main peak is pushed back to the ~ 7 MeV region, when the isoscalar 
mesons are included. Another information contained in Fig. ^| is that there is practically no 
contribution from the the spatial part — D p (ri, r^a^'cr 2 ) of the p-induced interaction. The 
dash-dotted curve (practically on the solid curve) corresponds to the calculation in which 
this part of the p exchange interaction has been omitted. 

The IVGDR strengths in 208 Pb, calculated with four different effective interactions, are 
shown in Fig.[5| In all four panels a strong collective shift from the free Hartree response 
to the full RRPA response is observed. This shift is predominantly caused be the isovector 
p meson field in the particle-hole channel. The isoscalar mesons pull the strengths slightly 
back to lower energies. The experimental IVGDR energy in 208 Pb is 13.5±0.2 MeV [31]. All 



the calculated peak energies are very close to this value, and it is difficult to discriminate 
among the different interactions. Although the IVGDR is mainly related to the coefficient 
of the symmetry energy, this relation is not very pronounced in the calculated peaks. For 
instance, the symmetry energy is 35 MeV for HS, and 43.7 MeV for NL1, but the peak and 
centroid energies of the corresponding strength distributions are very similar. 

In Fig.[| we display the transition densities for the IVGDR peak at 12.9 MeV in 208 Pb, 
calculated with the NL3 interaction. The proton and neutron contributions are plotted 
separately. The radial dependence is characteristic for the isovector giant dipole resonance: 
a node at the center, the proton and neutron densities oscillate with opposite phases in the 
surface region, they are in phase in the interior of the nucleus. 

Finally, in Fig.|5| we show the IVGDR strength distributions for four different spherical 
nuclei, calculated with the NL3 interaction. The main peaks are generally found very close 
to the experimental centroid energies, and the Landau width increases for lighter nuclei. 

D. Giant quadrupole resonances 

In a harmonic oscillator potential the microscopic structure of quadrupole excitations is 
simple: 2hu particle-hole configurations coupled by the residual interaction. In the isoscalar 
(isovector) channel the particle- hole interaction is attractive (repulsive), and it generates a 
collective state which is found in the energy region below (above) 2hu. In addition to the 
isoscalar giant quadrupole resonance (ISGQR) and isovector giant quadrupole resonance 
(IVGQR), part of the isoscalar strength can be found in the Ohu low-lying discrete states 
(seen subsection III. A). 

The ISGQR strength functions for four magic nuclei are shown in Fig|| The effective 
interaction is NL3. The collective downward shift from the Hartree response (dotted curves) 
to the RRPA response (solid lines), is clearly observed. The figures also illustrate the effect 



of the a/i-configurations (a empty state in the Dirac sea, h occupied state in the Fermi sea) 
on the RRPA response. As we have already explained, a rather large effect can be expected 
for isoscalar modes. Indeed, apart from 16 0, the positions of the main peaks are shifted to 
higher energy by 1-2 MeV by the inclusion of a/i-configurations. The effect, however, is not 
as large as in the monopole case. It also seems that the low-lying 2 + discrete states are less 
sensitive to the inclusion of ah configurations in the RRPA. 

The corresponding transition densities clearly exhibit the isoscalar nature of these exci- 
tations. In Fig.[7| we display the transition densities for the low-lying 2 + state, and for the 
isoscalar giant quadrupole resonance in 208 Pb: the neutron and proton densities oscillate in 
phase, the transition density of the discrete state is localized on the surface, while a large 
component in the interior region is observed for the transition density of the giant resonance. 

The isovector giant quadrupole resonance is less collective, because the repulsive isovector 
interaction mediated by the p meson is weaker than the attractive isoscalar interaction (er 
and uj meson fields). This leads to large Landau fragmentation of the strength distribution, 
similar to that observed in the IVGMR case. In Fig|| we display the IVGQR strength in 
208 Pb calculated with the NL3 parameter set of the effective Lagrangian. The distribution 
has a peak at 23 MeV, very close to the position of the experimental centroid energy at 
22 MeV |3"3fl . The transition densities of the state which corresponds to the peak energy 



are displayed in Fig. || As in the dipole case, the isovector nature of this mode is clearly 
observed in the surface region. 

E. Giant octupole resonances 

One of the best studied resonances is certainly the isoscalar giant octupole resonance 
(ISGOR) in heavy nuclei. In 208 Pb for instance, in addition to the low-lying 3~ state (cf. 
subsection III. A), the octupole strength distribution displays a pronounced fragmentation 
in two energy regions: the low-energy octupole resonance (LEOR) and the high-energy 
octupole resonance (HEOR) |32|]. The strength in the lower energy region corresponds 
mainly to Ihuo excitations, while Zhuj configurations dominate in the high-energy octupole 
resonance. The ISGOR strength distributions in 208 Pb are displayed in Fig.[H]. The effective 
interaction is NL3. In addition to the low- lying discrete state (represented as a vertical bar 
in arbitrary units), the LEOR strength is found between 6 and 10 MeV, and the HEOR 
strength is concentrated around 23 MeV. The RRPA results are similar to those obtained 
in non-relativistic Hartree-Fock-RPA calculations [§,[| , and compare well with experimental 



data 32 



The strength functions have been calculated with the response function method and, in 
addition to the full RRPA (solid curve), the Hartree response (dotted curve) is displayed, as 
well as the strength function calculated without the contribution of the ah configurations. 
Also for the isoscalar octupole mode a strong effect of the negative energy states is observed. 
Their contribution shifts the centroid of the HEOR to higher energy by few MeV. Without 
the inclusion of the ah configurations in the RRPA basis, the energy of the LEOR becomes 
even negative. This clearly shows the importance of working with a complete basis in the 
relativistic RPA framework. 
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IV. CONCLUSIONS 

Collective excitations in spherical nuclei have been investigated in the framework of a 
fully self-consistent relativistic random phase approximation (RRPA). By using effective 
Lagrangians which, in the mean-field approximation, provide an accurate description of 
ground-state properties, the strength distribution functions have been calculated both with 
the linear response method and by diagonalization of the RRPA matrices in configuration 
spaces. 

Results of RRPA calculations of multipole giant resonances and of low-lying collective 
states have been analyzed. For a set of spherical nuclei, the strength functions of normal 
parity multipole giant resonances: IVGMR, IVGDR, ISGQR, IVGQR, and ISGOR, have 
been calculated with the non-linear effective interactions NL1, NL3, NL-SH, TM1, and for 
comparison, also with the linear parameter set HS. The results have been compared with 
available experimental data. An excellent agreement has been found for the positions of 
the giant resonance peaks calculated with the non-linear effective interactions. The calcu- 
lated low-lying collective states are also in very good agreement with experimental data, 
particularly if one takes into account the fact that the effective interactions have not been 
in any way adjusted to properties of excited states. The relativistic mean-field Lagrangians 
with meson self-interaction terms provide not only an excellent description of ground state 
properties, but can be very successfully applied in the analysis of the dynamics of collective 
excitations within the RRPA framework. 

It has been shown that an RRPA calculation, consistent with the mean-field model in 
the no — sea approximation, necessitates configuration spaces that include both particle- hole 
pairs and pairs formed from occupied states and negative-energy states. The contributions 
from configurations built from occupied positive-energy states and negative-energy states 
are essential for current conservation and the decoupling of the spurious state. In addition, 
configurations which include negative-energy states have an especially pronounced effect 
on isoscalar excitation modes. For the isoscalar states the RRPA matrix contains large 
elements which couple the a/i-sector with the p/i-configurations. This is due to the fact 
that the contribution of the time-like component of the vector meson fields to these matrix 
elements, is strongly reduced with respect to the corresponding contribution of the isoscalar 
scalar meson field. The well known cancellation between the contributions of the a and 
u mesons, which for instance, takes place in the ground-state solution, does not occur in 
the RRPA matrix for the isoscalar states. The effect on the excitation energies of isoscalar 
modes can be up to few MeV. The isovector excitation modes, on the other hand, are 
much less affected by the contribution of negative energy states in isospin zero nuclei. The 
effect is slightly more pronounced in isospin non-zero nuclei. Isospin mixing in the isovector 
excitation modes has been observed in isospin non-zero nuclei, e.g. 208 Pb, and probably such 
effects could be rather important in studies of nuclei with a large neutron excess. 
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TABLES 
TABLE I. Parameter sets of the effective Lagrangians used in the present analysis. 





NL1 


NL3 


NL-SH 


TM1 


HS 


m [MeV] 


938.0 


939.0 


939.0 


938.0 


939.0 


TTv [Me V] 


492.25 


508.194 


526.0592 


511.198 


520.0 


m w [MeV] 


795.359 


782.501 


783.0 


783.0 


783.0 


777 p [MeV] 


763.0 


763.0 


763.0 


763.0 


770.0 


9a 


10.138 


10.217 


10.44355 


10.0289 


10.47 


9lu 


13.285 


12.868 


12.9451 


12.6139 


13.80 


9p 


4.975 


4.474 


4.3828 


4.6322 


4.04 


92 [far 1 ] 


-12.172 


-10.431 


-6.9099 


-7.2325 




53 


-36.265 


-28.885 


-15.8337 


0.6183 




C3 


0.0 


0.0 


0.0 


71.3075 




K nm [MeV] 


211.7 


271.8 


355.0 


281.0 


545.0 



TABLE II. Calculated and experimental excitation energies, and B(EL) values for the low-lying 
vibrational states in 208 Pb. The calculated values correspond to NL3 parameterization, the data 
are from Ref. [29|. Energies are in MeV, B(EL) values in e 2 fm 2i . 



L 77 

;r 

5" 

2+ 

4+ 



2.76 
3.26 
4.99 

4.95 



R 

2.61 

3.71 
4.07 
4.32 



B(EL) th 
499 xlO 3 
201 xlO 6 
2816 
998 xlO 4 



B(EL) e:r p 
(540±30)xl0 3 
330 xlO 6 
2965 
1287xl0 4 



TABLE III. The centroid energies of the IVGMR in ^ U8 Pb calculated in the RRPA with different 



effective interactions. The experimental centroid is from Ref. [32|. All energies are in MeV. 



NL3 



NL1 



NL-SH 



TM1 



Exp 



26.42 



26.96 



27.14 



26.28 



26.0±3.0 
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FIG. 1. The isovector monopole strength distribution in 208 Pb, calculated with the NL3 effec- 
tive interaction. 
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FIG. 2. The IVGDR strength function in 208 Pb, calculated with NL3. Solid curve: full RRPA 
response; dashed curve: the ah configurations are omitted from the RRPA basis; dash-double dot- 
ted curve: only the /^-mediated particle-hole interaction contributes to the RRPA matrix elements; 
dash-dotted curve: without the matrix elements of the spatial parts of vector meson fields. The 
free Hartree response (dotted curve) is displayed to illustrate the collective effect. 
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FIG. 3. The IVGDR strength distributions in 208 Pb, calculated with different parameter sets. 
Solid curve: full calculation; dash-double dotted curve: with only the p-mediated particle-hole 
interaction; dotted curve: free (Hartree) response. The arrows point to the experimental value of 
the IVGDR excitation energy. 
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FIG. 4. Transition densities for the isovector dipole state at 12.9 MeV in 208 Pb. The effective 
interaction is NL3. 
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FIG. 6. The ISGQR strength distributions in selected closed shell nuclei, calculated with NL3. 
Solid curve: full calculation; dotted curve: Hartree response; dashed curve: without the contribu- 
tion of the ah configurations. The vertical bars (in arbitrary units) denote the low-lying states. 
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FIG. 7. Transition densities for the low-lying 2 + state (a), and for the ISGQR (b), in 208 Pb 
calculated with the NL3 interaction. Dashed curves: proton densities, dotted curves: neutron 
densities, solid curves: total transition densities. 
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FIG. 8. The IVGQR strength distribution in 208 Pb calculated with NL3. The main peak is 
found at 23.0 MeV 
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FIG. 10. The ISGOR strength distributions in 208 Pb calculated with NL3. Solid curve: full 
calculation; dotted curve: Hartree response; dashed curve: without the contribution of the ah 
configurations. The vertical bar (in arbitrary units)denotes the position of the low-lying octupole 
state. 
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